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I. INTRODUCTION 
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In hadrons with a single heavy quark the dynamics of the hght degrees of freedom becomes independent of the 
heavy quark flavor and spin when the mass of the heavy quark is made arbitrarily large. This is known as heavy 
quark symmetry (HQS) SSl- This symmetry can be developed into an effective theory (HQET) 5] that allows 
a systematic, order by order, evaluation of corrections in inverse powers of the heavy quark mass. Ordinary HQS can 
not be applied directly to hadrons containing two heavy quarks. There, the kinetic energy term, needed to regulate 
infrared divergences, breaks heavy flavor symmetry [g!]. Only the spin symmetry for each of the heavy quark flavor is 
preserved. The symmetry that survives is heavy quark spin symmetry (HQSS), which amounts to the decoupling of 
the heavy quark spins for infinite heavy quark masses. In that limit one can consider the total spin of the two heavy 
quark subsystem (Sh) to be well defined. HQSS is sufficient to derive relations between form factors for the decay 
of hadrons containing two heavy quarks. That was first shown in Ref. JJ, where the authors adopted an approach 
where the two heavy quarks bind into a color anti-triplet which appears as apointlike color source to the light degrees 
of freedom. Applying the "superflavor" formalism of Georgi and Wise [1, allowed the matrix elements of the 

heavy-flavor-changing weak current to be evaluated between different baryon states. Semileptonic decays of the Be 
meson were also studied using HQSS in Ref. The formalism employed in [ll[ has been recently extended to 

describe semileptonic decays of be baryons to cc baryons [T3|. The scheme presented in does not rely on the 
"superflavor" formalism and HQSS is naturally implemented in it. In agreement with Ref. Q, the authors^ of Ref. [l^l 
found that spin symmetry for two heavy quarks enormously simplifies heavy to heavy semileptonic baryon transitions 
in the heavy quark limit and near the zero recoil point. As a result it is shown how an unique function, called the 
Isgur-Wise (IW) function, describes an entire family of decays involving doubly heavy baryons with total spin 1/2 and 
3/2. This imposes limitations to any quark model calculated form factors. Besides, the fact that all baryon matrix 
elements are given in term of just one function induces relations among different decay widths that, to our knowledge. 
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TABLE I: Quantum numbers of doubly heavy baryons analyzed in this study. J is the spin parity of the baryon, and Sh is 
the spin of the heavy degrees of freedom. I denotes a light u ox d quark. Mass predictions from Ref. ,14j] obtained using the 
ALl interquark potential of Ref.fisl] are also given. 



^ They find two differences with the results of Ref. Q, which cannot be eliminated by redefining the phases of the physical states. One 
difl^erence was already pointed out in [T3| . 
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have not been exploited before to check the vahdity of different quark model calculations. This is the main purpose 
of this letter. 

In a recent work [l^ we have studied, within a nonrelativistic quark model framework, static properties of doubly 
heavy baryons and their semilcptonic decays driven by the 6 — > c transition at the quark level. For the scmilcptonic 
decays we limited ourselves to spin 1/2 to spin 1/2 baryon transitions^. While we have shown that our wave functions 
have the correct limit for infinite heavy quark masses'^, we did not check HQSS constraints on the form factors or 
decay widths. Here we would like to extend our previous study on doubly heavy baryon b ^ c semilcptonic decays 
to include also doubly heavy spin 3/2 baryons and test our model and other s ag ainst HQSS predictions. These typ e 
of decays have been studied in different relativistic quark model approache s Il6l . [TtI [isj. with the use of HQET p^ . 
using QCD sum rules [l^l and three-point nonrelativistic QCD sum rules |2l|. or in the framework of the operator 
product expansion using the inverse heavy quark mass technique [2^ . Discrepancies between the results obtained 
in different quark model are sometimes very large. Therefore, it is worthwhile to use HQSS relations among decay 
widths to check the validity of the different calculations. 

In Table H] we summarize the quantum numbers of the doubly heavy baryons considered in this study"*. 



II. FORM FACTOR DECOMPOSITION 



Hadronic matrix elements can be parameterized in terms of form factors. For 1/2 — > 1/2 transitions the commonly 
used form factor decomposition reads 

(i?'(l/2), r' p' I F(0)7'^(l - 75)vl/^(0)| i3(l/2),r = uf:ip'){r (Fiiw) ~ -f^G^iw)) + (F^iw) - 75G2H) 

+v'^(Fs{w)-^5G,{w))}u^{p) (1) 

with \B{S)^r p) representing a baryon state with three-momentum p, total spin S, and spin third component r. 
The baryon states are normalized such that {B{S),r' p' \ B{S),r p) = (2'itY{Eb /ms) 5rr' 5^{p ~ p') being EB,mB 
the baryon energy and mass. The are dimensionless Dirac spinors, normalized as Ur'Ur — &„' ■ ^ v'^ are the 
four velocities of the initial and final baryon. The three vector Fi, F2, F^, and three axial Gi, G2, G3 form factors 
are functions of the velocity transfer to — v ■ v' or equivalently of the four momentum transfer {q = p — p') square 
q^ — m% + m^, — 2mBmB'UJ. In the decay [ ] ranges from uj = 1 [q^ = (Zmax = (™s ^ ) corresponding to 

zero recoil of the final baryon, to a maximum value given by w = Wmax = (™b + it^'b' ^ 'm^'i) / {'^mB'^B') \_q^ = n^f] , 
which depends on the transition, and where mi stands for the final charged lepton mass (we neglect neutrino masses). 

For 1/2 3/2 transitions we follow Llewellyn Smith [2^ to write 

{B'{3/2),r'p'\ ^\0)r{l~j,)^\0) \B{l/2),rp) = ufUp')T^^u^ (p) 

^ f^^g^.^ „ ,A^,^ ^ ^{g'^qp' - q'p'n + ^{g'^qp - <zV) + Q^(c.).9^'^) 75 
V rUB rrig nig J 

'Ga^H, 



{g'^i q'r) + ^{g'^qp' - q'p'n + C^{.)g'^ + ^q\A 
TUB rrig nig / 



(2) 



with p, p' the four-momenta of the initial, final baryon, and where we use the convention g^^ = (+, — , — , — ). u^^t is 
a dimensionless Rarita-Schwinger spinor normalized as u\ri = —5rr'- 

For 3/2 — > 1/2 transitions we have 

(i?'(l/2),r'p"| r(0)7''(l-75)*''(0)|i?(3/2),rp) = u?.' {p;')t>^ ^uf,{p) 

= f[T^^^{mB ^ mB',p ^p',q^ -q)]t/ (3) 



^ In that reference we missed a factor l/\/2 that affected our results for form factors. Decay widths were thus affected by a factor 2. An 
erratum has been sent. 

^ In the infinite heavy quark mass hmit the baryon should look like a meson composed of a light quark and a heavy diquark. 
* Note that the definitions of H;,;, and S^^ are interchanged in some references, with Hj.^. having Sh = and ^ having 5";, = 1. The same 
applies to and In tables we always quote the results corresponding to our convention (see Table 
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Finally for 3/2 3/2 transitfons, we believe there are 25 vector plus 25 axial form factors. The amount of form 
factors suggest a different strategy in this case and thus we do not show the form factor decomposition. 

In Ref . jl4l| we presented results for 1/2 1/2 transition form factors. In a similar way one can evaluate all form 
factors for 1/2 < — > 3/2 transitions. It is not the purpose of this work to present results for all individual form factors. 
Instead, we would like to study to what extend they obey the restrictions imposed by HQSS. 



III. HQSS CONSTRAINTS ON FORM FACTORS FOR SEMILEPTONIC DOUBLY HEAVY BARYON 

DECAY 

We quote in what follows the results obtained in Ref.[l^, using HQSS and near zero recoil, for the semileptonic 
be — *■ cc baryon decay with the initial baryon at rest. There it was found that all hadronic matrix elements were given 
in terms of just one universal function (77(0;)), known as 

1 
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71 

where here B stands for a S or baryons. The IW function which controls the S decays is different to that appearing 
in Q decays since the IW function depends on the light degrees of freedom. The IW function 77 is approximately one 
at zero recoil {ri{LU = 1) w 1), as can be deduced from vector conservation in the limit of degenerate b and c quarks^. 

Similar results can be obtained for semileptonic bb — > be baryon decays, but there will have different IW functions 
because of the heavy flavor symmetry (HFS) breaking in hadrons with two heavy quarks. 

Let us see the implications of the above relations for the form factors calculated in quark models. 
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the IW function. Indeed HQSS predicts 
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A. 1/2 —» 1/2 transitions 

Near zero recoil the three vector structures 7'^, v'^ and u'^ present in Eq.([T]) give, up to corrections proportional 
to \q \ that cancel near zero recoil, the same contribution. On the other hand, the Dirac's structure of the axial form 
factors G2 and G3, and due to the anti-diagonal nature of 75, give contributions that are again proportional to |g | 
and thus cancel near zero recoil. To the extent that for the actual heavy quark masses we are close enough to the 
infinite heavy quark mass limit, Eqs. (|4l5p imply the following restrictions on form factors 

Bfcc-Bcc -^(Fi+F2 + F3) = ^Gi=77 (10) 

B'bc^B,, {Fi+F2+F3)=0 ; -^Gi^r, (11) 

The same relations can be derived respectively for B^b — Bbc a-nd B^b B'f^^ decays. 

In Fig [1] we show our results for the above quantities evaluated for the Sf,c — > Sec, SJ^^ —f 'E.cc and Sbf, Sf,c, 



^bb ^'bc transitions. We have used the ALl interquark potential of Ref. jl5j and actual heavy quark masses (baryon 



masses are given in Table |T|. We do not show very similar results obtained for transitions involving fl baryons. We 



^ Note the authors of Ref. tl2| missed a global normalization factor I/V2. An erratum has been sent. 




see in the figure that the above restrictions are, to a good approximation, satisfied by our calculation over the entire 
Lo region. 



B. 1/2 < > 3/2 transitions 



In this case one can see that all contributions generated by the ^ CY , , Cq , and Cq form factors are propor- 
tional to |g I, cancelling thus near zero recoil. The same happens for the 9^7'' dependence of the form factor and 
the q'^p''^ dependence of the form factor. On the other hand the g'^^ dependence of the axial part of the current 
survives near zero recoil. The restrictions imposed by Eqs. (|6l7l8p are in this case 



Bbc - i?:, - \/l ( - "^^c-c ) + i^B,^ Eb.. - ) + cA^V (12) 

BL - Sc*c - ^ [-^^"'^L " ^B'J + -^{^biEb:^ - ml.J + C^^ = ^ (13) 



BL - Sec - V —i^B',^ - "^Sc. ) - ~ ruB^EBj + = ^ (14) 

For Bbb B^^, B^f^ — > Bt,c and B'^^^ the relations obtained are given respectively, and with obvious changes, 

by Egs. p^ . and again ^1)1 but in the latter case with the factor changed to 1/^/2. 

We show now In Figl2]our results for the S^c S*^, SJ,^ — > E!*^ and S^^ S*^ transitions, and for the 'Rhh S^^, 
SJjj — Sbc and S^^^ — > S[,^ transitions, again evaluated with the ALl interquark potential of Ref . [l5l| and actual heavy 
quark masses. We do not show results for U, baryons which are very similar to the ones presented. We see again our 
calculation is in accordance with HQSS constraints. 



C. 3/2 —> 3/2 transitions 



In this case we have not evaluated explicitly individual form factors. Here we proceed as follows: we evaluate the 
hadronic matrix elements in Eq. ([9]) for different spin configurations, selecting only vector or axial components for 
which the matrix element does not cancel near zero recoil. The corresponding matrix elements are also evaluated in 
the quark model of Ref. 1J|. By comparison of the two calculations we get the IW function. In this way one can 




FIG. 2: Left panel: relations in Eqs. (|12ll3ll4p for Etc — >■ S*^, Hf,^ — * H*c and H^c ^ Hcc transitions. We also show, for 
better comparison the results already shown in Fig. [1] Right panel: similar relations for E,bb —> H^^., Etc and 31^ 5';,^, 

transitions. All the results have been obtained using the ALl interquark potential of Ref . [l5l| . 




FIG. 3: Left panel: Different IW functions obtained for Hj^^, — > H^c transitions (black curves) using the vector or the axial part 
of the weak transition current, and for different spin configurations. For better comparison we also show the corresponding 

results obtained for 1/2 — > 1/2 and 1/2 < > 3/2 transitions. All results have been obtained using the ALl interquark potential 

of Ref. [Tsl . Right panel: same as left panel for bb — > be transitions. 



obtain five different functions, two of them with the vector part of the current and three others with the axial part. In 
the infinite heavy quark mass hmit these five functions should coincide among themselves and with the ones obtained 
in 1/2 ^ 1/2 and 1/2 < — > 3/2 transitions. 

The results are shown in Fig. [31 What we see is, that to a good approximation, better in the bb be case as one is 
closer to the infinite heavy quark mass limit, all 1/2 — > 1/2, 1/2 < — > 3/2 and 3/2 3/2 transitions are governed in 
terms of just one function. As mention before this function is different for the be ~> cc and bb ~* be cases due to HFS 
breaking. 



IV. SEMILEPTONIC DECAY 



In this section we present our results for semileptonic decay widths and compare them with the ones obtained 
in other quark model approaches. In some cases there are large discrepancies between different calculations. The 
fact that, at least for the bb be case, we are not far from the infinite heavy quark mass limit suggests that some 
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calculations might be inconsistent with HQSS. 
The decay width is given by 



du! \/u 



1 /:^"Wu. 



(15) 



where Gp = 1.16637(1) x 10 MeV ^[2J| is the Fermi decay constant and \Vch\ is the modulus of the corresponding 
Cabibbo-Kobayashi-Maskawa matrix element. is the leptonic tensor defined as 



(fk (Pk 



2E 2E 



- 5^^\q -k-k') {k'^'k'' + k'^kf" - g^"'k ■ k' + ie^""^'^k'^kp) 



(16) 



where fc, k' represent the momenta of the final charged lepton and antincutrino respectively. We use the convention 



^0123 _ 



1. Using Lorentz covariance one can write 



2^ 9^9'^ 



(17) 



where neglecting neutrino masses 



2e 



with 



(18) 



(19) 



Note that for a light lepton / = e,/i we can neglect terms in rn^/q^ over most of the q^ interval and thus use 
B{q^) « -A{q^). 

The hadron tensor is given by 



x{b\S'),t' p' vI/^(0)7.(/-75)^''(0) B{S),rp 



(20) 

In Ref. T3| it is shown how 1/2 1/2 hadronic matrix elements are evaluated within our model. The extension to 
the 1/2 < — > 3/2 and 3/2 3/2 cases is straightforward. 

In Table HIl we compare our results for S — s- S transitions with the ones calculated in different models. Our central 
values have been obtained with the ALl potential of Ref. [15] , while the errors shown indicate the spread of the results 
when using four other interquark potentials, three more taken from Ref. [T5| and another one from Ref. psj . In all 
cases we have used a value |V"c&| — 0.0413. Our results are in a global reasonable agreement with the ones in Ref. [l^ 
where they use a relativistic quark model evaluated in the quark-diquark approximation®. In Ref. p^ . and using 
HQET, results around a factor of 4 larger than ours are found^. We believe this factor of 4 discrepancy stems from the 
fact that in Ref. ^ISi] the author approximates 'q{ijj) by ri{l) in the calculation of the decay width*. In Ref. [2l| they 
obtain results similar to the ones in the previous reference^. Finally in Ref. [l3| they obtain in general much larger 



Note the results we show under Ref. Il^ are a factor of 2 smaller than the originally published. The reason beint that the authors of 
that reference also missed a normalization factor l/\/2 for diquarks with two equal quarks |2^ . An erratum has been sent. 
Note the results we show under Ref. jl9l | are a factor of 2 larger than the originally published. There is a factor \/2 wrong in the 
normalization of matrix elements that affects the published results j27| . 

If we take for instance the approximate expression in the Eq. II21I I , which is closer to the approximations used in Ref. jl9| , and make 
ri{u}) = 77(1) we get r(Hi,c — > Sec) = 9.4 X 10~^* GeV in agreement with the result in Ref. |19| . On the other hand if we take the actual 
ri{uj) values we get r(Hi,(, S^c) = 2.4 X lO"^* GeV, roughly a factor of 4 smaller and in agreement with our full calculation result. 
We must say we believe the calculation in Ref. |21| is affected by the same normalization mistake made in the original calculation in 
Ref. [T9I I as they give Fi + F2 + = rj instead of Fi + F2 + F3 = y/2T]. To our understanding their present results have to be multiplied 
by a factor of 2. 
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TABLE II: Decay widths in units of 10"^'* GeV for doubly heavy H baryon semileptonic decay. Our central results have been 
obtained with the ALl potential of Ref. [isl ]. The errors show the spread of results when using four other interquark potentials 
taken from Refs. [isl. [2a|. We have used a value! Vet | = 0.0413. I stands for a light charged lepton, I = e, fi. For resuhs with ^ 
and * see text for details. The results of Ref. [2j| are given as quoted in Ref. |16i |. 
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TABLE III: Same as Table [IT] for doubly heavy Q baryon semileptonic decay. Decay widths are given in units of 10"^* GeV. 



results for all transitions. In this latter calculation the authors take the Bethe-Salpeter equation model to analyze 
the weak transition matrix elements between two heavy diquarks, and then use "superflavor" symmetry [1, 
to evaluate the transition matrix elements at the baryon level. The global results show a contradiction between the 
calculation by Guo et al. (TtI in one hand, and ours and the one by Ebert et al. [l^ on the other. 

In Tabic [TTTl we show results for ^ transitions. Again we get a global reasonable agreement with the calculation 
by Ebert et al. [ll|. 

It is worthwhile to mention that for the case of baryons with a he heavy quark content the actual physical states 
S and will be an admixture of S^c, Sj,^ and fibc, ^'i,c respectively. If we look for instance at our model predictions 



we see the widths are very different for transitions involving or 
decay widths could thus give information on the admixtures. 



and fJftc or fi^^^. Accurate measurements of 



V. HQSS CONSTRAINTS ON SEMILEPTONIC DECAY WIDTHS 

To the extent that one is close enough to the infinite heavy quark mass limit and near zero recoil we can combine 
the HQSS results in Eqs.(|4][9l) with Eq. fTTl) . to approximate the tensor product £^'''71^^ by 
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and similar ones for bb — > be decays. 

Working in the strict near zero recoil approximation, oj « 1 or equivalently q^ quite close to its maximum value 
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and Aiff') w —B{q^) near g^^x- these circumstances, and using 

HQSS predicts that the different decay widths are in the relative ratios 
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4:3: 2:3: 1:5 (30) 

In Table HVl we show the above ratios obtained in different models. Our results and the ones by Ebert et al.fl^ are 
in reasonable agreement with the HQSS predictions in this strict near zero recoil approximation. On the other hand 
the results by Guo et al.\v^ deviate heavily form the above predictions. This disagreement does not improve much 
by using a different decay width to normalize the ratios 
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TABLE IV: Decay width ratios for semileptonic be — » cc and bb — > be decay of doubly heavy H and baryons compared to 
the HQSS predictions in the strict near zero recoil approximation. Our results have been obtained with the ALl potential of 
Ref. [151] . / stands for a light charged lepton, I = e, fj.. 

We can relax the strict near near zero recoil approximation to obtain more accurate predictions based on HQSS in 
the following way. For the actual doubly heavy baryon masses LOmax ~ 1-22 (1.08) for be — > cc (66 — > 6c) transitions 
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while the different differential decay widths dT/duj show a maximum at around oj « 1.05 (1.01). We can thus still use 
w « 1 and A{q'^) « —B{q^). On the other hand the quantities {v ■ jq^ ■, (v' ■ qY /q^, {v ■ q){v' ■ q)/q^-, that are all 
equal to 1 near zero recoil, can deviate rapidly from 1 because of the q^ factor in the denominator. What is true, in 
and around the maximum of the differential decay width, is that we can reasonable approximate 

{v-qf W ■q){vq) 



q- 
jv' ■ qf 
g2 



{v' ■q){v- q) 



q- 



(31) 



With the above consideration we can still predict approximate ratios between different decay widths that one expects 
to be satisfied to an accuracy of 20 ^ 30%. We have chosen to define those ratios so that they are all equal to one, 

(32) 
(33) 
(34) 
(35) 



nBL - 




. nsib- 




w 1 


3r(Bfc,- 


-^B*jDi) ^ 


^ 3r(B*,- 


Bbc Ivi) 




^B*JPi) ^ 


T{Bbb- 




« 1 




Bcclvi) 


1 r(i?66 




r(Sfce - 


^ Bcclvi) 


, nsib- 


-> Bbc Ivi) 


w 1 


ir(BL- 


Bcclyi) 


\ ^{Bbb 


-^B^ln) 


nsi^^Btjvi) 






nB*bb 


B*bM) 


T{Bbc ^ B^Jui) + ^T{Bt,c ' 


-^B*Jui) 


TiBbb ^ 


Bbclvi) + 


\ V{Bbb ^ Bl^ iDi) 



Note, we consider as independent the phase-space integrals of ri'^{uj)A{q'^) and ?7^(w)A((7^)i^^-^^^. 

In Tables IVIVII we show the above ratios evaluated in different quark model approaches. Once again calculations in 
this work and the ones in Ref. |.16i] are compatible, within the expected accuracy, with the approximate ratios obtained 
using HQSS results. On the other hand the deviations found in the results by Guo et al. are, in most cases, too large. 





Tliis work 
E Q 


[16] 
H Q. 


[17] 

H n 


nB'bc- 

iV{Bbc- 


^B*M) 
-^B*Jvi) 




1.041°;°? 1.04-0.03 


0.79 0.82 


0.68 — 


T{Bbc- 

'inB'bc 


-BlJvi) 
-^Bcclvi) 




0.82l°;°? 0.84l°;« 


1.22 1.17 


2.72 — 


nBi- 
\nBL- 


^Bcclvi) 
-^Bccl^l) 




0.94+°" 0.97l°;j; 


1.28 1.26 


10.6 — 


nBL- 


-BlJvi) 




0.89+°'" 0.941°; J? 


1.01 1.01 


1.08 — 


T{Bbc^Bccln)^ 


2 ^{Bbc~ 


>B*jDi) 



TABLE V; Decay width ratios for semileptonic be — > cc decay of doubly heavy H and Q baryons. In all cases the approximate 
result obtained using HQSS is 1. Our central results have been obtained with the ALl potential of Ref. [l5| . The errors show 
the spread of results when using four other interquark potentials taken from Refs. [isl. l25l]. I stands for a light charged lepton, 
I — e, 



VI. SUMMARY 



We have checked the constraints imposed by HQSS on form factors and decay widths. To our knowledge those 
constraints have not been exploited before to check the consistency of different quark model calculations. We have 
shown that our calculation is consistent with HQSS. The ratios in Eqs. (|29|30p . obtained using HQSS with strict zero 
recoil approximation, and the approximate ratios in Eqs. (|32j|35[). w here we have relaxed that approximation, compare 
well with the results in our model and the one by Ebert et al. [l6[ , but they are incompatible with the calculation in 
Ref. ■ We think that although this is not enough guarantee for the predictions here and in Ref. [l^ to be fully 
correct (in fact the few results in Refs. [l^, [2l| are not incompatible with HQSS constraints while they are a factor of 
four larger than ours), it certainly indicates problems either in the model or in the calculation performed in Ref. p^j . 
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This work 


fl6] 


[17j 

H n 




^B'bM) 
-^Bbclvi) 






0.99 


0.99 


U.Uo — 


r( Bu, — 


*^bc '■^l) 
--BlWi) 




0.86+o"5 


0.96 


0.99 


y.Oo — 


nB*bb- 

3 r(Shb- 


^Bbclvi) 
-^B',Jvi) 


0.98«»^ 




1.01 


1.03 


36.4 — 


m*bb- 


^B*bJPi) 




0.93^°:1J 


1.01 


1.01 


0.31 — 


T{Bbb^Bbc m)+\ T{Bbb^B;^ Wi) 



TABLE VI: Same as Table |V] for semileptonic bb — > be decay of doubly heavy H and Q baryons. 
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